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Abstract 

The forward-backward asymmetry in top pair production at Tevatron has been reconfirmed by 
the CDF collaboration with 5.3 fb~^ of accumulated data. These measurements also report that 
the asymmetry is the largest in regions of high invariant mass Mjf and rapidity difference |Ay|. We 
consider light colored sextet scalars appearing in a particular non-supersymmetric 50 (10) grand 
unification model within the 126 scalar representation. These scalar states have masses in the range 
of 300 GeV — 2 TeV consistent with the requirements of gauge coupling unification and bounds on 
the proton lifetime. The cross section and the total asymmetry can be simultaneously explained 
with the contributions of these scalars within la. We find that the simultaneous fitting of the cross 
section, the total asymmetry and the asymmetries in different rapidity and M^i bins gives only a 
marginal improvement over the SM contribution. We also study various production mechanisms 
of these colored sextet scalars at the LHC. 
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I. INTRODUCTION 



The CDF collaboration had measured the forward-backward asymmetry, ApB, in pair 
production in the ti rest frame in 2008 at Tevatron with 3.2 fb~^ of collected data, as [l| 

A^fcos ^ > 0) - iVfcos e < 0) , , , , 

Afb = -T77 ^ ^77 ^ 7 = 0.193 ± 0.065 stat ± 0.024 syst 1 

iV(cos e > 0) + iV(cos ^ < 0) ^ ' ' ^' 

where 6 is the scattering angle of top quark in the ti rest frame. This result had been 
confirmed by DO collaboration based on 0.9 fb~^ of integrated luminosity They reported 
Afb = 0.19±0.09(stat)±0.02(syst) and Afb = 0.12±0.08(stat) ±0.01(syst) for exclusive 4- 
jet and inclusive 4-jet events respectively consistent with CDF results. These measurements 
have attracted a lot of attention due to more than 2a deviation from the Standard Model 
(SM) predicted value of = 0.058 ± 0.009 Q. In the SM, the Afb identically vanishes 
at leading order (LO). However, at next to leading order (NLO) in QCD, it can arise from 
(a) the interference between tree level SM amplitude and the box diagram, (b) radiative 
corrections to qq annihilation and (c) interference between different amplitudes contributing 
to gluon-quark scattering. Several independent New Physics (NP) scenarios have been 
advanced to explain this discrepancy. 

Recently, CDF have presented new results using 5.3 fb~^ of data sets in which Afb is 
reported to be 0.158 ± 0.074 (stat+syst). Along with this new value of Afb, they find 
interesting dependences of Afb on the invariant mass of the ti pair and their rapidity 
dependence. The asymmetry is more prominent in the large invariant mass region of > 
450 GeV with more than 3a deviation and in the large rapidity difference |Ay| > 1 region 



with around 2a deviation from the SM predicted value [10|. On the other hand, some 
other observables related to ti pair production at Tevatron show good agreement with the 
SM predicated values. The measured parton level ti -pair production cross section a^^^ = 
agrees with the SM predicted value of cr|*'^ = 7A5toll calculated with MCFM 



7.70±0.52 
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Similarly the experimentally measured invariant mass distribution is also consistent 



with the prediction of the SM at NLO pj[. Hence, while discussing the new physics scenarios 
to explain the new results, namely the mass and rapidity dependence of Afb, they must 
not introduce large corrections to either the total ti cross section au or the invariant mass 
distribution Mu- Some recent attempts 0, [isj have been made in order to explain these 
new observables along with the updated measurements of Afb and au- 

The standard model gauge structure {SU{3)c x SU{2)l x U{1)y) allows a finite number 
of different representations of scalar particles which can couple an up or down quark to the 
top quark. The possible cases include a set of colored octet, singlet, triplet and sextet scalars 
each for uu ti and dd ti processes. 

(8,2,i), (l,2,i), (3,3,|), (6,3,|) ioiuu/dd^ti 

(3,1, 1), (6,1, 1) for MM ^tt only (2) 

(3, l,i), (6,l,i) for ciJ -> only 
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The effects of these scalars on ApB have been studied in a model independent way in several 
papers. For example, the results of Shu et al show that colored sextet and triplet scalars 



are able to explain the anomaly while the analysis of Jung et al [6| favours the singlets and 
sextets. On the other hand, the results of Arhrib et al [tI] show that sextet diquarks could 
not fit the Aps and cross section simultaneously within la. However all these analyses were 
based on old observations oi ApB- According to the new CDF data, the central value oi ApB 
has significantly come down. Also the distributional preferences of ApB in invariant mass 
and rapidity have been reported. So it is an interesting exercise to study these scalars in the 
light of new observations. In the present study we investigate the light colored sextet scalars 
appearing in a well motivated 5*0(10) grand unified theories as a possible explanation of 
ApB well as new observables simultaneously based on current data. 

Colored scalar fields naturally emerge in a well motivated class of grand unified the- 
ories. For example, the representations (8,2,|) and (3, 1,|) reside in a 45-dimensional 
Higgs field of SU{b). It is interesting to note that in any simple renormalizable version of 
non-supersymmetric S't/(5), the 45 Higgs together with 5-dimensional Higgs is necessarily 
required to generate viable masses of charged fermions |l6|. It has been shown through 
detailed studies in reference js], ^ that both the scalar states (8, 2, |) and (3, 1, |) can have 
masses in the range of 300 GeV - 1 TeV consistent with the requirements of gauge coupling 
unification and bound on the proton lifetime. In addition, the contribution of colored triplet 
scalar to the production of ti at the Tevatron can enhance the forward-backward asymme- 
try and account for the experimental result without spoiling the successful standard model 
prediction for the total cross section. We investigate a similar possibilities in more predic- 
tive and attractive class of grand unified theories based on the S'O(IO) gauge group. The 
remarkable feature of 50(10) is that its 16-dimensional irreducible spinor representation 
accommodates a complete family of fermions, including the right-handed neutrino. This 
complete unification of quarks and leptons opens up the possibly of connections between the 
charge fermions and the neutrino sector. Furthermore, 5*0(10) has the left-right symmetry 
group SU(2)l X SU{2)ji as a subgroup, making the implementation of the both the type-I 
and the type-II seesaw mechanisms very natural in these theories. The 45 dimensional scalar 
representation of SU{5) resides in both the 126 and the 120 dimensional scalar representa- 
tions of 5*0(10) which can couple with ordinary fermions through the Yukawa interactions. 



The 126 Higgs field plays an important role in gauge symmetry breaking [17| as well as it 



is essential for viable fermion masses [18|]. However the color triplets of 126 couple with 
fermions through symmetric leptoquark couplings and can induce the rapid proton decay 
if assumed light. On the other hand, the 120 Higgs field has antisymmetric coupling with 
fermions but it is not required if one sticks to the minimal Higgs content of the theory. 
We show in this work that the 126 Higgs has diquark colored sextets and octets at TeV 
scale consistent with gauge coupling unification and proton decay bounds and study the 
role of the sextet scalars as the possible candidates to explain the anomaly in ti production 
observables. 

This paper is organized as follows. In the next section, we study the scalar spectrum of a 
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particular SO{10) model and discuss the constraints coming from gauge coupling unification 
and proton decay. In section III, we study the role of light colored sextet scalars on the tt 
pair production observables. In section IV, we will study the signatures of these scalars at 
both 7 TeV and 14 TeV LHC Finally, we summarize our results in section V. 



II. LIGHT COLORED SCALARS IN 50(10) MODEL 

We consider non-supersymmetric 5*0(10) as a basic framework of our model. It has 
been pointed out in recent studies jl3] that an adjoint 45-dimensional scalar representation 
(x) of 5*0(10) together with one 16 or 126 (S) Higgs can govern the entire breaking of 
50(10) gauge symmetry down to the SM. If one sticks to the renormalizable version of the 
seesaw mechanisms then the representation 126 is indispensable, since it breaks the SU{2)r 
group and gives neutrino masses through seesaw mechanisms. In addition, one needs 10 



dimensional Higgs {(f)) to obtain a realistic fermion mass spectrum [18[. We have given the 
decompositions and full SM spectrum of these scalar fields in Table IIVI in Appendix. Note 
that X contains two SM singlets (xs, Xs) and S contains one (Sg) SM singlet that acquire 
vevs at GUT scale and break 5*0(10) to the SM group. The 10 and 126 Higgs contain the 
SM doublets (02^ ^2) and (S2, S2) respectively which can mix through a renormalizable term 
XijXki^ijkim(pm in the scalar potential. For consistent fermion mass spectrum, one has to 
keep (at least) one linear combination of these doublets light upto the electroweak scale, 
which plays the role of SM Higgs doublet and triggers the electroweak symmetry breaking. 
This requires a fine tuning in the parameters of the Higgs potential. Assuming such fine 
tuning in parameters, a detailed numerical analysis has been carried out recently for viable 



fermion mass spectrum in this model in Ref. p^. It has been shown that such model 
can provide very predictive structure of fermion masses if a global U{1)pq (Peccei-Quinn) 
symmetry is imposed and can produce realistic fermion mass spectrum which is in excellent 
agreement with the present data extrapolated at the GUT scale. 

We consider a non-supersymmetric 50(10) framework with the minimal Higgs fields 
10 + 45 + 126 in our attempt to explain the forward-backward asymmetry in tt production 
at Tevatron. Following previous studies, we assume that only one linear combination of the 
weak doublets of 10 and 126 remains light and becomes the SM Higgs. Further, we also need 
to assume that the scalar submultiplets which can potentially contribute to the asymmetry 
in the production of top quarks also remain light in the range of 300 GeV - 2 TeV. Among the 
possible options allowed by the SM gauge symmetry shown in Eq. |2l the 126 contains three 
sextets 23(6, 1, |), E4(6, 1, |), Si2(6, 3, — |) and a pair of a octets Si5(8, 2, |), Si5(8, 2, — |). 
These fields couple to the 16-plet matter through Yukawa interactions. Furthermore, all 
these fields are diquark (have only quark-quark coupling) in nature and do not mediate the 
proton decay. Some other components of the 126 scalar (like Si,S7 and S13) also have 
the correct quantum numbers to influence ti production but they have leptoquark coupling 
which induce rapid proton decay if assumed light. If no artificial suppression via Yukawa 
couplings is arranged their masses should not be below than 10^^ GeV due to proton decay 
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constraints. 

The interaction of 126 Higgs field to the 16- dimensional matter fields can be written 
in its most general form as 
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-Cy = ^FijijjBC ^TpTgTrTsTti^jEpgrst (3) 

where the indices i,j denote family indices, p,q,.. = 1,..,10 are 5*0(10) indices, C is the 
Dirac charge conjugation matrix and B = TiTsT^T^Tg is the equivalent of the charge conju- 
gation matrix for 5*0(10) that ensures the invariance under 5*0(10). Fj's are representations 
of the Clifford algebra associated with the Lie algebra of the 5*0(10) group and are given 
. F is the Yukawa coupling matrix and it is symmetric by its 5*0(10) properties. 
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After the decomposition |2lf| of Eq. [H the couplings of S3, S4, S12 and S15 to matter can 
be written as 

V2F., {dl{T\bC-'d%t.li + ul{TXbC-'d%f.~,,^) (4) 

where = ^A"^ and {A = 1,..,8) are the Gell-Mann matrices of 5f/(3). a,b,c are 
color indices. Clearly, all these fields have the right couplings to influence the asymmetry 
we are interested in. However, in order to be relevant for asymmetry at Tevatron, all these 
fields or at least one of them must be sufficiently light. On the other hand, such light 
states will contribute to the running of gauge couplings and hence viability of their being 
light is constrained by the unification of gauge couplings and present bound on the proton 
lifetime. We thus show that it is possible to achieve light colored scalars with successful 
gauge coupling unification in a consistent way in our model. 

In the absence of any new particle thresholds between the weak and GUT scales, the 
running of gauge couplings at one-loop level is given by 



-1 -1.,. X bi f Mgut\ 



(5) 



where acvT represents the gauge coupling at the unification scale Mqut- are the ap- 



propriate one-loop /3 function coefficients [2^ and 2 = 1,2,3 stands for f/(l)y, SU{2)l and 
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5*?7(3)c respectively. Their values for the SM with one light Higgs doublet are hi — 
62 = and 63 = —7. It is easy to check that these values for hi do not unify the gauge 
couplings since SM does not predict gauge coupling unification in the first place. The pres- 
ence of new particles between weak and GUT scale can change the running and it can be 
easily incorporated by replacing hi in Eq. [5] with effective one-loop coefficients Bi defined 



by [23 



B^ = ,^^Y.b^^MMcEim (6) 
V HMgut/Mz) ^ ' 
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where bu is the one-loop coefficient of the additional particle / of mass Mj lying between 
Mz and Mqut- Following Giveon et al. {23I, Eq- [5] with contributions from Eq. [6] can 
provide successful gauge coupling unification at one loop level if they satisfy following two 
conditions: 

'^^"'^^-"/"^= 0.716 ±0.005 (7) 



-B23 


_ B2 


-B; 


B12 


Bi 


-B2 


B12 


= B, 


-B2 



8 3/8-sin^^v^ 
167r 3/8 - sin^ 184.9 ± 0.2 



5a hi{MGUT/Mz) HMgut/Mz) 



(8) 



We use the present experimental measurements of the SM parameters 2J] to derive the 
above numbers. In any given model, Bij depend only on the particle content and associated 
mass spectrum and conditions ([7]) allow us to constrain the mass spectrum of particles that 
leads to an exact unification at GUT scale. We give a list of the different submultiplets of 
10, 126 and 45 scalar representations and the corresponding contributions to coefficients Bij 
in Table fllVp in the Appendix. 

In order to present a consistent analysis, we now discuss the constraints coming from 
proton decay. In nonsupersymmetric GUTs, this process is mediated by baryon number 
violating gauge interactions which induce a set of effective dimension six operators at low 
energies that conserve B — L. In the 50(10) scenario we consider here, such gauge bosons 
are integrated out at the GUT scale {mx,Y = Mqut) and therefore proton decay constrains 
Mgut from below. The most stringent bounds coming from the latest experimental limit 



on partial decay lifetime of proton Tp {p — )■ vr^e"'") > 8.2 x 10'^'^ years 25|] implies 



Mgut ^ (mja^^rrp)^ > 2.3 x 10^^ GeV, (9) 

where rup = 0.938 GeV is the proton mass. Some of the submultiplets of 10 and 126 Higgs 
are leptoquark scalars (for example, 0i,Si,S6 and etc.) and are associated with d = 6 
proton decay operators. We suppress their contribution to proton decay by making them 
super heavy ~ Mgut as we will explain in next paragraph. 

As mentioned earlier in this section, it is necessary that at least one of the submultiplets 
(S3, S4, S12 and S15) of 126 remains light in order to explain the forward-backward asym- 
metry in top quark production. Typically in theories with two or more widely different mass 
scales, if a submultiplet of a full Higgs multiplet acquires a vev ~ M, the members of that 



multiplet acquire a mass ~ M [261] . Any scenario which differs from this would require some 
fine tunings in the parameters of the scalar potential. To check the viability of such fine 
tunings, a complete detailed analysis of scalar potential minimization and its diagonalization 
is required. However this is beyond the scope of present work and we assume that such fine 
tuning is possible in our case. In order to avoid further unnecessary fine tunings, we assume 
that the remaining submultiplets of the scalar fields 10,45,126 are super heavy and have 
natural masses of order Mgut- In other words, we assume that only those submultiplets of 
scalar fields remain light and have masses Mj in between the weak scale and the GUT scale 
and may potentially contribute to the forward-backward asymmetry in tt production. 
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With all these considerations, we now check the compatibility of light colored scalar 
states of our interest with the unification of the gauge couplings and constraints on GUT 
scale coming from proton decay bounds. Following the strategy of j^, we determine an upper 
bound on GUT scale at the one loop level assuming that any one of E3, E4, E12 and E15 is 
responsible for asymmetry and is accordingly in the mass range of 300 GeV - 2 TeV. For this, 
we numerically maximize Mqut while imposing the condition that the solution satisfies Eq. 
[3 The additional constraints we put on the solution are 300 GeV < rriY,^, 171-^^,171^:^^2,171^-^^ = 
rrij^-^^ < Mgut and Mqut < Mpianck = 10^^ GeV. The results of our numerical analysis are 
shown in Fig. [1] and Fig. [2l We get viable gauge coupling unification consistent with proton 
decay limits in two different scenarios. 




500 1000 1500 2000 



FIG. 1. The maximum value of Mgut obtained for different values of ms^j by assuming S3 light. 
The dashed line stands for the lower bound on Mgut due to the proton lifetime. Viable gauge 
coupling unification is achieved in the region between ?tt-5]i2=300 GeV and the dashed line. 

(A) We get successful gauge coupling unification for the scalar diquark E3 having mass in the 
range of 300 GeV to 2 TeV as shown in Fig. [H The light E3 also requires light Ei2(6, 3, — |) 
and there is a clear correlation between their masses. There exists an upper bound on ms^j 
for a given value of m^g. For example, when = 600 GeV we have m^ij ^ 2 TeV. The 
other two scalar states remain heavy, namely, E15 ~ 10^ — 10^^ GeV and E4 ~ Mgut- 

(B) Unification of gauge coupling is also achieved with light sextet diquark state E4(6, 1, |) 
as shown in Fig. [2j Unlike light E3 in the previous case it does not require any other light 
submultiplet at TeV scale. The maximum value of the GUT scale does not change apprecia- 
bly with and stays well above the present proton decay limits shown by the dashed line 
in Fig. O Successful unification in this case requires only one state at intermediate scale 
^^212 ~ 10*^ GeV and E3, E15 and remain superheavy (~ Mgut)- 

From the results of the detailed analysis carried out in this section, we conclude that either 
E3 or E4 can remain light and influence the forward-backward asymmetry in ti production 
through the processes uu — )■ ti and dd — )■ ti respectively. We do not get viable gauge coupling 
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FIG. 2. The maximum value of Mqut obtained by assuming S4 light. The dashed line stands for 
the lower bound on Mqut due to the proton lifetime. Viable gauge coupling unification is achieved 
in the region between two lines. 

unification with a light colored octet state S15 and hence this case will be ignored in further 
analysis. Note that we have presented consistent unification analysis at the one loop level. 
The allowed masses of S3 and S4 would change slightly if one considers two-loop effects in 
the running of gauge couplings. However they would still remain within the TeV range. 

III. COLORED SEXTETS AND FORWARD- BACKWARD ASYMMETRY OF 
TOP QUARKS 

The light colored sextet scalars E3 and S4 contribute to ti pair production through u- 
channel exchange as shown in Fig. [3l E3 interferes with the SM contributions for uu and 
cc initial parton states while S4 interferes with the SM contributions for dd , ss and bb initial 
parton states. The contributions of initial parton states cc, ss and bb to the overall process 
pp — )■ ti will be suppressed due to their small parton distribution functions (PDF). However 
we include all these contributions in our analysis. 

Let us denote incoming quark momentum by Pg, incoming anti-quark momentum by Pq, 
outgoing top momentum by pt and outgoing anti-top momentum by pt with the following 
definitions : 

a/I 

p,,5=^(l,0,0,±l), (10) 
Pt,t = [Et,±\^t\{sme,0, cose)], (11) 
= Et[l,±Pt{sm6,0,cos9)] (12) 

where Et = l3t = \pt\/Et = a/1 — 4m^ / s and 9 is the angle between momenta of the 

incoming quark and the outgoing top quark in the center of mass (cm) frame of the partons. 
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FIG. 3. Contributions from light sextet scalars to the tt production at the Tevatron. 



Neglecting masses of all quarks except the top, the Mandelstam variables in parton cm are 
defined as follows : 



S = {Pq+Pqf = {pt+PtY 

i = {Pq - Ptf = {Pq - Ptf 



X1X2S 
2 



Acos6'), 

u = (Pq- Ptf = (Pq - Ptf = mi-'-{l + (3t cos 9) 



(13) 
(14) 

(15) 



where s is the cm energy of proton and antiproton in laboratory frame, Xi and X2 are the 
fractions of momentum carried by the partons inside proton and antiproton respectively. 

With these notations and conventions, the matrix amplitude squared (averaged and 
summed over initial and final color and spin indices respectively) for [qq — )■ ti ) can be 
written as follows : 

,2 2sl 



total 
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l + Msin2^ + cos2e 



\fi3 



u,d\2 



9 u — m% 



I fU,d\i 
Iil3 I 



[l + PtCosef + 



Pt cos ( 



2 , 



(16) 



where /"g and /{'g are related with the original coupling F13 of Eq. lUby the following relation 



(17) 



'13- 

For our numerical study, we have used the leading order PDF sets of CTEQ6L [27|J to 
convolute with the partonic cross section to obtain hadronic cross section. We set our 
renormalization and factorization scale to fiR = fip = rrit. The top mass is taken to be 
rrit = 172.5 GeV at which we also evaluate strong coupling = 0.1085. We use i^— factor 
of 1.3 to rescale our LO results for a{tt) to match with NLO QCD prediction 28 . 

We calculate the total cross section a(ti ), ApB as defined in Eq. [H as well as ApB in 
lAFl > 1, lAFl < 1, Mtt < 450 GeV and Mu > 450 GeV, where lAFl is the difference 
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of top and anti-top quark rapidities i.e., |Ay| = Yt — Y[ in tt rest frame. The present 
experimentally measured values of all these observables, their values predicted in the SM 
and corresponding contributions needed from NP are listed in Table [B 



Observables 


Experimentally 
Measured Values 


SM Contribution 


Contribution needed 
from NP 


Cross section 
Afb 

Afb{Mu > 450GeV) 
Afb{Mu < 450GeV) 
AfbH^YI > 1) 
Afb{\AY\ < 1) 


7.70 ± 0.52 
0.158 ±0.074 
0.475 ±0.112 
-0.116 ±0.153 
0.611 ±0.256 
0.026 ±0.118 


0.058 ± 0.009 
0.088 ±0.0013 
0.04 ± 0.006 
0.123 ±0.018 
0.039 ± 0.006 


0.1 ±0.083 
0.387 ±0.1133 
-0.156 ±0.159 
0.488 ± 0.274 
-0.013 ±0.124 



TABLE I. The observables with their experimentally measured values, their values predicted in 
the SM and corresponding contributions needed from NP. The contributions needed from NP are 
obtained by subtracting the SM contributions from experimentally measured values. 

We perform a analysis to simultaneously fit all the observables shown in Table HI For 
this, we define the following function 

i=l ^ ^ 

where the sum runs over all the six observable quantities. Pj's are the theoretically calculated 
values of these quantities as a function of couplings and masses of scalars in our model and 
Oj's are the mean values of these observables. (Ti's denote la errors in Oj. The is 
numerically minimized to obtain the best fit over all six observables. The more robust 
statistic to quantify the quality of fit is reduced-x^ which is defined as the x"^/^ where u is 
number of degrees of freedom (d.o.f.) in the analysis. For the SM, the value of total is 
17.26 and the value of xV^^ is 2.88. 

We now present a detailed numerical analysis for the contributions of S3 and S4 sepa- 
rately. 

A. Diquark (6, 1, |) 

In Fig. m we plot the cross section and the forward backward asymmetry for tt production 
at Tevatron as a function of the coupling /"g for four different masses of the colored sextet 
scalar E3. In showing the contribution from new physics, we subtract the SM contribution 
from the experimentally measured value of Afb- 

From Fig. HJ we see that sextet S3 of mass 300 GeV can barely satisfy both a** and Afb 
constraints for a very narrow range of coupling /"g that too at 2a. For large masses S3 can 
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FIG. 4. The tt production cross section (left) and the forward backward asymmetry (right) as a 
function of the couphng /fg for masses msg = 300 GeV (Red), 900 GeV (Green), 1500 GeV (Blue) 
and 2100 GeV (Black). The unshaded and the shaded region correspond to la and 2a experimental 
bounds respectively. 





mY.; = SOOGeV 


jTiEa = 900GeV 


= 1.5TeV 


msg = 2.1TeV 


Observables 


Fit pull 


Fit pull 


Fit pull 


Fit pull 


Cross section 
Afb 

Afb{Mu > 450GeV) 
Afb{Mu < 450GeV) 
Afb{\AY\ > 1) 
AfbHAYI < 1) 


8.0404 0.6546 
0.0237 -0.919 
0.0336 -3.119 
0.0453 -1.616 
0.0154 1.078 
0.0161 0.234 


8.2761 1.1078 
0.0617 -0.462 
0.0968 -2.5614 
0.126 -1.3196 
0.0282 1.1585 
0.0358 0.3932 


8.3004 1.1546 
0.0690 -0.3729 
0.1113 -2.434 
0.1446 -1.253 
0.0285 1.1602 
0.038 0.414 


8.3076 1.1685 
0.0715 -0.3429 
0.1163 -2.389 
0.151 -1.23 
0.0284 1.159 
0.0391 0.420 




14.83 


11.24 


10.48 


10.22 




3.71 


2.81 


2.62 


2.55 


\fn\ 


0.549 


1.319 


2.105 


2.905 



TABLE II. Results of analysis carried out for different values of mxig . The best fitted values for 
each observables along with their respective pulls are shown. The pull measures the deviation in 
the fitted value of the observable from its mean value. For NP contributions, the number of d.o.f. 
is 4 (No. of observables — No. of parameters.) 

satisfy both the constraints for large range of couphng /^g within lo" of the experimental 
bound. 

The results of the analysis are shown in Table [Tll We show the best-fit values of 
all the observables along with their respective pulls. The minimum values of and the 
corresponding values of parameter /"g obtained at the minimum are shown for different 
masses of Eg. The overall fits get better with increase in m^g. For all masses, we get the 
largest pulls corresponding to Afb in the Mti > 450 GeV region where it gives more than 
2a deviation. All the other observables can be fitted within 1.2a. Although the total for 
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FIG. 5. The tt invariant mass distribution for NLO-SM, and for various masses of colored sextet 
S3 for the best fitted values from our analysis. 



NP contribution is better than for the SM, the x^/^ values are worse than the SM value 
for smaller values of S3 mass. The x^/ shows slight improvement relative to the SM only 
for masses greater than 1.5 TeV. Hence, the sextet S3 in our model can satisfy the total 
cross section and the total asymmetry within la while it is incompatible with asymmetries 
in the large invariant mass region and the large rapidity region for the same parameter space. 



An empirical relation between f^^ and m^g is obtained and can be written in approximate 
form as : 

l/^l =0.148 + 1.31^^. (19) 

I J 131 ^ rp^y V J 



Another important constraint in the ti production comes from the invariant mass distribu- 
tion of ti pair. This distribution has been measured by CDF collaboration and is shown in 
Fig. |5]for various values of S3 masses with CDF data and SM-NLO prediction. We use the 
best fit values of coupling /"g for various masses as shown in Table HTl for evaluating the con- 
tribution of NP to the Mtt distribution. The SM contribution to the da/dMu distribution 
in Fig. [5] has been evaluated to the full NLO order as given in Ref. 0]. While evaluating 
the contributions of NP and its interference with the SM to invariant Mu distribution, we 
multiply the contribution with K-factor of 1.3. However, it is highly desirable to include 
full NLO corrections to NP to make more reliable prediction on the invariant mass distri- 
bution. We see that lower values of S3 masses fit the distribution better than larger values 
of S3 masses. 
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FIG. 6. The tt production cross section (left) and the forward backward asymmetry (right) as a 
function of the couphng ff^ for masses m^^ = 300 GeV (Red), 900 GeV (Green), 1500 GeV (Blue) 
and 2100 GeV (Black). The unshaded and the shaded region correspond to la and 2a experimental 
bounds respectively. 

B. Diquark (6, 1, i) 

Next, we study the diquark S4 to look at its effect on ti pair production at the Tevatron. 
The contribution of E4 to tt production has been shown in Fig. [3] and proceeds through a 
dd initial state. In Fig. El we plot the cross section and the forward backward asymmetry 
for tt production at Tevatron as a function of coupling /fg for four different masses of colored 
sextet scalar S4. As stated earlier, in showing the contribution from new physics, we subtract 
the SM contribution from the experimentally measured value of ApB- 

Because of the fact that the PDF for the d-quark is smaller than that for the u-quark, we 
need larger values of the coupling ff^ to generate the contribution to the observables of tt 
pair production. From Fig. El just like for S3, S4 of mass 300 GeV can barely satisfy both 
cr** and Afb for a very narrow range of coupling ff^ that too at 2cr but with different range 
of coupling. For larger masses S4 can satisfy both the constraints for large range of coupling 
ff^ within lo" of experimental bound but with a wider and different range of couplings as 
compared to S3. 

The results of the analysis for S4 are shown in Table UTTl Similar to the previous case, 
the overall fits get better with increase in rrij^^. As seen earlier, the x'^/'^ values are worse 
than the SM value for smaller values of S4 mass. Similar to the previous case, the best-fit 
relation between ff^ and can be put into approximate form as : 

If^J = 0.273 + 2.48^^. (20) 

The perturbativity argument regarding the strength of a generic coupling requires 
g^/4:7r < 1 , which allows, in principle, the coupling to be as large as ~ 3.5. In Figs. HI 
and El we show the cross section and the Afb up to /"g'^ = 4. However, these are not the 
canonical couplings which enter in the Lagrangian of Eq. |H /"g*^ is related to the canonical 
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"^E4 = 


300GeV 


"^E4 = 


900GeV 


"^E4 = 


1.5TeV 


rriY.^ = 


2.1TeV 


Observables 


Fit 


pull 


Fit 


pull 


Fit 


pull 


Fit 


pull 


Cross section 


8.0157 


0.6071 


8.2649 


1.0863 


8.2734 


1.1027 


8.2968 


1.1478 


Afb 


0.0230 


-0.9277 


0.0604 


-0.4775 


0.0490 


-0.6143 


0.0710 


-0.3499 


Afb{Mu > 450GeV) 


0.0299 


-3.1515 


0.0985 


-2.546 


0.1041 


-2.4972 


0.1192 


-2.3636 


Afb{Mu < 450GeY) 


0.0428 


-1.6249 


0.1299 


-1.3070 


0.1365 


-1.2827 


0.1547 


-1.2164 


Afb{\^Y\ > 1) 


0.0184 


1.0971 


0.0345 


1.201 


0.0338 


1.1934 


0.0350 


1.2012 


Afb{\AY\ < 1) 


0.0173 


0.2441 


0.0405 


0.4313 


0.0405 


0.4313 


0.0429 


0.4513 


x' 


15.07 


11.23 


11.09 


10.15 




3.78 


2.81 


2.78 


2.54 




1.041 


2.508 


3.875 


5.527 



TABLE III. Results of analysis carried out for different values of ms^ . The best fitted values for 
each observables along with their respective pulls are shown. The pull measures the deviation in 
the fitted value of the observable from its mean value, v denote number of degree of freedom. For 
NP contributions, the number of degree of freedom is 4 (No. of observables — No. of parameters. 



coupling Fi3 according to the Eq. [T71 Using this relation, we find that the values of the 
couplings fif which we use in our analysis satisfy the perturbativity. 

The invariant mass distribution of ti pair corresponding to contribution of E4 is shown 
in Fig. [7] for various values of S3 masses with CDF data and SM-NLO prediction. We use 
the best fit values of coupling ff^ for various masses as shown in Table UTT] for evaluating the 
contribution of NP to the Mti distribution. From the fig., we see that all values of S3 masses 
are more compatible with the distribution in the large Mu bin and fit the distribution better 
than S3. However, there is a little tension in the distribution for bins 450 GeV-500 GeV 
and 550 GeV-600 GeV. 

We now discuss the constraints on masses and couplings of the colored sextet S3 and S4. 
These constraints have been discussed in detail in Ref. jsol where authors have analyzed 
the electroweak precision data (EWPD) to obtain lower bound on the mass of the sextet 
scalars. They find that EWPD does not give a lower bound much above 100 GeV. The 
constraint is weak because there is no custodial SU{2)c violation. The most robust bound 
on sextet masses comes from direct search of these scalars at LEP-11 putting a lower limit 
of 105 GeV on their masses. At Tevatron, the most stringent bound comes SjJ] from the 
search of narrow resonances in the dijet mass spectrum. They reported lower mass bound 
for diquark to be 290 GeV. The S3 can produce same sign dileptons through decay into 
two top quarks while S4 safely avoids same sign dilepton constraints. The constraints on 
the mass of S3 from the search of same sign dilepton signature is however weaker than the 
bound which comes from the search of narrow resonances. 

The other stringent constraints come from low energy processes such as — mixing. 
The contributions of the sextet colored scalars to this mixing has been studied in some detail 
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FIG. 7. The tt invariant mass distribution for NLO-SM, and for various masses of colored sextet 
S4 for the best fitted values from our analysis. 



in Ref. 
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E3 contributes in — mixing at tree level while S4 contributes through a 
box diagram. The bound on parameter for the S3 is |-Re(/;^2/iT)l ^ 5.76 x 10~^ for itih = 1 
TeV. The bound on coupling of S4 is lA^/nP < 1-7 x 10"^° for nin = 1 TeV. However, 
these bounds can be relaxed if the couplings to second generation is minimized which we do 
in our analysis for search these scalars at LHC. 



IV. SEXTET DIQUARKS S3 AND S4 AT LHC 



35 



The detailed phenomenology of sextet diquarks has already been performed in Refs. 33 



where they found that such scalars can be discovered at the LHC with masses around 
few GeV to 2 TeV. We have already shown in previous sections that S3 and S4 can have 
masses around this mass range in order to achieve unification at the GUT scale and explain 
the anomaly in the ti production forward backward asymmetry at Tevatron. 

The colored sextet scalar diquarks can be produced in the following channels at LHC: 

1. Resonant production in s-channel : pp— )■ S3/S4 t + u/d + X , 

2. Pair-production : pp-^ S3S*/S4S* + X ^ {t + u/d) + (t + u/d) + X, 

3. Single production with top quarks : pp— t- T^'i^^i/T,*^ + X. 



The channel 1 has been explored in detail in Refs. [33|,|3J|. The matrix amplitude squared 
and the resulting parton level cross section have been given in these references. The full 
NLO QCD corrections to qq annihilation to scalar sextet diquark resonant state has been 
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performed in Ref. [3^. In this channel, uu and ud initial states dominate over uu and ud 
initial states because of the large PDFs of quarks compared to antiquarks. Also, it is inter- 
esting to compare the cross sections corresponding to S3 and S4 production. The S3 and 
S4 get contributions from uu and ud initial state respectively. The S3 gets enhancement 
due to large u-quark PDFs while S4 gets enhancement from two sources : (a) due to com- 
binatorics from initial state, the luminosity of S4 is du ® ud while that of S3 is uu, and (b) 
from the relations [11] and [201 it can be seen that the coupling ff^ is larger than /^g. Hence, 
the cross section for S4 is almost 5 times larger than S3 production cross section for low 
masses and is about 1.5 times larger for large values of S's masses. The best strategy to 
discover S3 4 in this channel would be to determine the invariant mass distribution of t+j 
and look for narrow resonances of S3 4 as discussed in detail in Ref. 
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The channel 2 has been explored in detail in Ref. 35| . The production process is mediated 
through QCD interactions through gg fusion and qq annihilation and hence depends only on 
sextet masses. The matrix amplitude squared and the resulting parton level cross section 
have been given in this reference. The channel in which S3 decays to tttt has been analyzed 
in great detail in Ref. 35| for 14 TeV LHC. They propose a reconstruction in the multijet 
plus same-sign dilepton with missing transverse energy samples to search for tttt final states 
from sextet scalar production. The decays of S4 would yield {t + j) + [t + j) which can be 
probed in 8-j channel of which two are b-jets and all jets are hard jets. The cross sections 
for both S3 and S4 are large enough so that they can be discovered in lower mass range at 
14 TeV LHC. In Fig. [HI we show production cross sections for channel 1 and 2 for various 
cm energies of LHC and for various possible initial states. To calculate cross sections, we 
evaluate couplings /"3 and ff^ from relations [19] and [20] respectively and assume f^^ = f^i 
and ft, = ft,. 




0.01 



0.001 



400 600 800 1000 1200 1400 1600 1800 2000 
Mass of sextet (in GeV) 



10000 


1 


1 1 


1 1 1 1 1 1 

^3 ■ 

Vs=7TeV ^4 - 


1000 






uu+cc^ ~- — - 


100 








10 






Dii+cc 


1 








0.1 


1 




1 1 1 1 1 1 


400 


600 800 


1000 1200 1400 1600 1800 2000 



Mass of sextet (in GeV) 
5^3,45^3 4^ + X (left) and s-channel resonant 



FIG. 8. The cross section for pair production pp 
production (right) of pp — )• 113^4 — )• t + u/ d + X with various initial states at the LHC. The values of 
the couplings f,^ and ff^^ are evaluated using relation [19] and [20] respectively. We assume couplings 
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The S3 and S4 sextets can also be produced in association with antitop quarks. The 
cross sections for production of Ss+t and T^^+i pair have been shown in Fig. [9l Because 
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FIG. 9. The cross section for process pp 
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of the large couplings of S4, the cross section for S4 production is comparable to that for 
S3 production. We see that for the interesting mass range of the sextets, the cross section 
is of the order of SM ti cross section which makes this channel very promising. The search 
strategy in this channel would be to look for tij final states and search for resonances in the 
invariant mass of light jet with antitop quark. 



V. SUMMARY 

In this paper, we have investigated the role of colored scalars as a possible explanation of 
large forward-backward asymmetry in ti pair production at Tevatron. We consider a partic- 
ular non-supersymmetric 5*0(10) model where such scalars reside in 126 dimensional scalar 
representation which plays a crucial part in GUT symmetry breaking and also generates 
small neutrino masses through seesaw mechanism. We find that either S3 or S4 colored 
sextet submultiplet of 126 can remain light and provide viable gauge coupling unification 
consistent with the present bounds on proton decay. 

Colored scalars in the context of forward-backward asymmetry at Tevatron have been 
studied in Refs. All these papers are based on old CDF data and do not include 

new observations. In this paper, we show that the contributions of light colored sextet 
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scalars S3 and E4 of mass > 300 GeV can enhance the forward-backward asymmetry with- 
out spoihng the a(tt) and the invariant mass distribution. We perform a analysis to 
simultaneously fit all the new observables along with total asymmetry and total cross sec- 
tion and find the best fitted values of the coupling for various masses. We find that all the 
observables can be fitted within 1.3a of experiment errors except for the Aps ^tt > 450 
GeV region where we can fit the observation at 2.3a. From the analysis, we conclude that 
the colored sextet scalars of masses 1.5 TeV-2.1 TeV can provide a marginal improvement 
over the SM observations if all observations are simultaneously considered in the fit. 

In this paper, we focus on the study of light colored sextet scalars in 50(10) model and 
their effects to the yli?^ in the light of new observations reported by CDF. We have shown 
that such scalars can emerge in particular S'O(IO) model having masses of the order of TeV 
scale. However the detailed analysis carried out in context of new and old CDF data are 
applicable to any sextet scalar. 

We also discuss the various production mechanisms of S3 and S4 at the LHC and find 
that these scalars will have observable cross section to be discovered in future. These scalars 
can be produced in pairs owing to pure QCD interactions and this channels is promising 
at LHC rather than Tevatron. They can also be produced in s-channel resonance and then 
can be probed in their decay to t + j events. Also, the other promising channel is to search 
them in single production in association with anti top quarks. The best strategy to search 
for these scalars in all these channels would be to look for narrow resonance in invariant 
mass of top quarks with light jets. 
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Note Added: While we were finishing the present work, Ref. 15| was appeared dealing with 



same kind of colored sextet state (6, 1, —3). Our results for S3 scalar qualitatively agree 
with their results. Compared to that work, we also study another sextet state and show 
that such scalars naturally emerge from theoretically well motivated 5*0(10) model. 



VI. APPENDIX 

In this appendix, we list all the sub-multiplets of the 10(0), 126(S) and 45(x) dimen- 
sional scalar representations of 50(10) and their contributions to Bij = Bi — Bj coefficients 
(where Bi, i = 1,2,3 are one loop (3 function coefficients for f/(l)y, 5f/(2)i and 5f/(3)c 
respectively). We also present these sub-multiplets in terms of their Pati-Salam subgroup 
(5?7(4) X SU{2)l X SU{2)ji) notations. {R—,RO,R+) represents components of the field 
which is triplet under SU{2)ji. The indices of the doublet of 5f/(2)i(5f/(2)ij) are denoted 
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by = l,2(a',/3' = l',2'). The index of the fundamental 4-plet of S'?7(4) is denoted by 
fi{iy) = 4 where fi.i' = 1, 2, 3 represents the SU{?>) subgroup indices. 



Fields {SU{2,),SU{2),Y) 


Pati-Salam Notations 


AB23 


AB12 


0i(3,l,-i), 0i(3,l,i) 


0M4, 0^' 


1 
6 


1 

15 


02(1,2, i), 02(1, 2,-i) 


0al', 02' 


1 
6 


1 

15 


XI (8, 1,0) 


x^ 


1 
2 





X2(3,l,i), X2(3,l,-i) 


X/2 , X4~ 


1 

6 


1 

15 


X3,8(l,l,0) 


^(15)^ ^(M) 








X4(3,2,-|), X4(3,2,|) 


XAi4a2', X^i^ol' 


1 
6 


1 

3 


X5(3,2,i), ;e5(3,2,-i) 


XflAal', XflPa2' 


1 
6 


7 

15 


X6(l,3,0) 


XaP 


1 
3 


1 

3 


X7(l,l,l), X7(l,l,-1) 







1 
5 


5^i(3, 1, — |), Si(3, 1, |) 




1 
6 


1 

15 


S2(l,2,l), S2(l,2,-i) 




1 
6 


1 

15 


^3(6,1, 1) 
S4(6,l,l) 
S5(6,l,-|) 
^6(3, 1, |) 
5^7(3,1,-1) 
5^8(3,1,-1) 


yim 

y{R-) 
^pv 

y{R+) 
^pi 

y{RO) 
^pi 

yiR-) 
^P4 

y{R+) 

^44 


5 
6 
5 
6 
5 
6 
1 
6 
1 
6 
1 
6 


32 
15 
2 

15 
8 

15 
4 

15 

1 

15 
1 

15 


^9(1, 1,0) 








5^10(1,1,-1) 


y{RO) 
^44 





1 

15 


5^11(1,1,-2) 


yiR-) 
^44 





1 

15 


Si2(6,3,-1) 


^pu 

^ ap 


3 
2 


18 
5 


Si3(3, 3, |) 


yP4: 
^ al3 


3 
2 


9 
5 


Si4(l,3,l) 


V44 

^ a/3 


2 
3 


1 

15 


5^15(8, 2, ^), 1115(8, 2, -^) 


y P V 

'^/i q1'' "^P a2' 


2 

3 


8 
15 


5^16(3,2,1), Si6(3,2,-1) 


y 4 V 

a2" ^4 al' 


1 
6 


7 

15 


5^17(3, 2, |), Si7(3, 2, -|) 


y 4 V 

^p al" ^4 a2' 


1 
6 


17 
15 



TABLE IV. Different sub-multiplets of 10(0), 126(S) and 45(x) dimensional scalar representations 
of 50(10) and their contribution to Bij coefficients of Eq. [H Various notations used are explained 
in text. 
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